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Are small-scale sub-structures a universal property of 
galaxy halos? The case of the giant elliptical NGC 5128 
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ABSTRACT 

We present an analysis of the spatial and chemical sub-structures in a remote halo 
held in the nearby giant elliptical galaxy Centaurus A (NGC 5128), situated ~ 38 kpc 
from the centre of the galaxy. The observations were taken with the Advanced Camera 
for Surveys instrument on board the Hubble Space Telescope, and reach down to the 
horizontal branch. In this relatively small 3.8 kpc x 3.8 kpc field, after correcting for 
Poisson noise, we do not find any statistically strong evidence for the presence of small- 
scale sub-structures in the stellar spatial distribution on scales > 100 pc. However, we 
do detect the presence of significant small spatial-scale inhomogeneitics in the stellar 
median metallicity over the surveyed field. We argue that these localized chemical 
substructures could be associated with not-fully mixed debris from the disruption of 
low mass systems. NGC 5128 joins the ranks of the late-type spiral galaxies the Milky 
Way, for which the stellar halo appears to be dominated by small-scale spatial sub- 
structures, and NGC 891, where localized metallicity variations have been detected 
in the inner extra-planar regions. This suggests that the presence of small-scale sub- 
structures may be a generic property of stellar halos of large galaxies. 
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"-H . 1 INTRODUCTION 



The early history of galaxy assembly leaves its signatures 
in the properties of the stellar content of the outer regions 
of galaxies. The currently-favored scenario of galaxy forma- 
tion in a hierarchical context predicts that the vast major- 
ity of stars in the outskirts of galaxies should be accreted 
from disrupted sate llite galaxies (e.g. lAbadi et al. 1 120061 : 
IZolotov et aTl 12009). The ages, metallicities, spatial distri- 
bution, and amount of sub-structures in the faint outskirts 
of present-day galaxies are therefore directly related to their 
assembly histories, and to issues such as the suppression of 
star formation in small halos. 

Recent wide field observations of the Milky Way have 
led to the identification of a number of large-scale stel- 
lar features, attribut ed either to past accretion of iden- 
tified satellites (e.g. Ilbata et all Il994l: lYannv etafl |200(1 
l2003l ; iNewberg etafl 120021 : iMaiewski et all |2002| ), others 



1 This work was based on observations with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 



of unknown progenit ors (e.g. iGrillmair fc Dionatos 1 l200d : 
iBelokurov et al.l 120071). and others to the disruption of glob - 
ular clusters (e.g. lOdenkirchen et a l. 2001; Grillmair 1 2009'). 
The accretion of stars from satellites is clearly a contribu- 
tor to the stellar haloes of galaxies. Quantitative analysis of 
wide-field imaging data from the Sloan Digital Sky Survey 
indicates that the stellar halo of the Milk y Way is highly 
structured on small scales (|Bell et al.ll2008l ). 

For the nearest Milky Way analogue NGC 891, the anal- 
ysis of deep and high resolution imaging data sampling a 
quadrant running parallel to the high surface brightness disk 
of the galaxy, and extending outward to more than 10 kpc 
from the plane of the galaxy, identified strong evidence for 
the presence of highly signific ant small-scale sub structures in 
the median metallicity map jlbata et all 120091 ). The degree 
of substructures present on the stellar density distribution 
in the halo of NGC 891 is found to b e strikingly comp arable 
to what is present in the Milky Way ijlbata et al.ll2009l ) . The 
extra-planar regions of NGC 891 appear to be composed of 
a large number of incompletely-mixed sub-populations. The 
discovery of a giant stream that loops around the galaxy 
and a series of low surface brightness features in its out- 
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skirts d emonstrates that NGC 891 had an active accretion 
history (|Mouhcine et alJlioiOah . 

Being the nearest giant elliptical galaxy, NGC 5128 
and its outer regions in particular have been targeted ex- 
tensively as a unique testing ground for stellar population 
and formation models of elliptical galaxies JSoria et al ]| 19961; 
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Harris et all Il999l; [Marleau et al.~l l200d; Harris fc Harris I 
200d . 120021 : llleikuba et al. Il200ll l2002l l2005h . The colour- 
magnitude diagrams of the stellar populations in the remote 
regions of the galaxy suggest that those regions are domi- 
nated by moderately metal-rich stars, i.e., [M/H] ~ —0.65, 
with a very wide red giant branch, suggesting that stars at 
large radii span a wide rang e of metallicities (|Harris et al.l 
Il99d ; IHarris fe Harris 1 120021 ). Different stellar age indica- 
tors, i.e., the red clump, the asymptotic giant branch bump, 
the period distribution of long period variables, globular 
clusters, suggest that the outer regions of NGC 5128 are pop- 
ulated by old stars, with a likely presence of a small fraction 
of intermediate age stars (Soria et al.| [l996; Rci kuba et al. I 
l200ll , l2003l ; IWoodlev et alj|20ich . 



Possible formation histories of elliptical galaxies come 
generically in t hree flavours, f i rst an early and ra pid in-situ 
formation (e.g. lLarson 1 1 1974 IHarris et allll995l ). secondly 
later mergers of pre-existin g disk galaxies (e.g. iToomrel 
Il977l ; I Ashman fc Zepf Hl992T ). and finally as the results of 
multiple dissipationless merg ers and accretion s of distinct 
protogalactic fragments (e.g. ICote et al.l 119981 ) . Significant 
effort has been devoted over the years to test these scenar- 
ios using a variety of discriminants, e.g. stellar ages, chem- 
ical abundances, kinematics, and globular cluster systems. 
These different formation histories should give rise also to 
different stellar halo structures. The in-situ formation would 
predict relatively few substructures, as the formation epoch 
was many dynamical times ago. The stellar halos of elliptical 
galaxies that have grown hierarchically through mergers and 
accretions would be however substantially lumpy, and pos- 
sibly as structured as those predicted by the most extreme 
cold-dark matter hierarchical formation models. 

The purpose of the present paper is to examine the dis- 
tribution of the stellar populations in the outer regions of 
NGC 5128, and to investigate in particular the dumpiness 
of the spatial and chemical structure of the halo population 
of a giant early-type galaxy for the first time. The data upon 
which this study is b ased have been presented in detail in 
iReikuba et al. (2005, hereafter referred to as R05). The lay- 
out of this paper is as follows: in ij2]we briefly introduce the 
data set used in the paper, while JJ] presents the analysis to 
uncover small-scale sub-structures over the surveyed area. 
The implications of our findings are discussed in §4\ and we 
draw our conclusions in 

Throughout this paper we assume a dist ance of 3.8 Mpc, 
i.e., distance modulus o f (m — M)o = 27.92 l|Reikuba et al. I 
l2005l ; IHarris et al.ll2009h . We note that NGC 5128 is located 
at relatively low Galactic latitude (£ = 309.52°, b = 19.42°), 
and therefore the targeted field suffers from significant 
(though not large) extinction from foreground dust. Follow- 
ing R05, we adopt a reddening value of E(B — V) = 0.11 
i|Burstein fc Heileslll982l ; ISchlegel et al.lll998h . which corre- 
sponds approximately to Ay — 0.36 and Ai = 0.21. 
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Figure 1. The location of the ACS field studied in this contri- 
bution lies ~ 38 kpc to the South of NGC 5128, superimposed on 
a photographic image extracted from the Digitized Sky Survey. 
The scale of the region is shown both in arcmin and kpc at the 
distance of N GC 5128. The X and Y axes of the ACS mosaic, as 
constructed bv lReikuba et al~l J2005T) are displayed on the image; 
these correspond approximately to the directions North and East, 
respectively. 



2 DATA 

As described in detail in R05, a single deep pointing was 
observed with the Advanced Camera for Surveys (ACS) 
camera on-board the Hubble Space Telescope (HST) in 
the halo of the galaxy NGC 5128 (program GO-9373), at 
a location approximately ~ 38 kpc South of the nucleus 
(Fig. [IJ. The target field was chosen to avoid any pe- 
culiaritie s of the galaxy, su c h as jet-induced star- forming 
regions dMould et al. I l200d ; IReikuba et al~l I2002T), shells 



l|Malin et al.lll983T l. and dust lanes (|Stickel et al.ll2004h . The 
F606W and F814W passbands were used (approximately 
broad-band Johnson V and Kron-Cousins /), with a total ex- 
posure of 8.58 hr (12 full-orbits) in each filter, whic h reaches 
F606 W~ 30, F814W- 29. ALLFRAME software (|Stetson I 
1 19941 ) has been used to measure the photometry of all de- 
tected sources on the ACS images. To select stellar objects, 
we applied the criteria of R05 (summarised in their Fig. 3), 
namely those sources with ALLFRAME sharpness parame- 
ter within their defined hyperbolic envelope, and with chi- 
squared < 3. The left panel of figure [2] shows the foreground 
reddening-corrected colour-magnitude diagram (CMD) of all 
objects detected in our images and classified as stars in the 
deep ACS halo field. The CMD are shown as density con- 
tours to reveal features in otherwise crowded regions. In ad- 
dition to the wide sequence of red giant stars, red clump 
stars are visible at I a ~ 27.5 and (V — J) D ~ 1.0, and those 
on the AGB bump at J D ~ 26.5 and (V - I) ~ 1.3. 

Careful analysis of the completeness corrections and 
photometric uncertainties has been conducted via Monte- 
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Figure 2. The left panel shows the foreground reddening-corrected colour-magnitude diagram of the targeted halo field. The dashed line 
shows the 50% detection completeness level. The right panel shows the colour-magnitude diagram of foreground stars as predicted by 
the Besangon Galactic population model in the direction of NGC 5128 over the same field-of-view as covered with a single ACS pointing. 
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Figure 3. Panel 'a' shows the surface density of RGB stars in the ACS field in the coordinate system of Fig. 1, binned into a 25 X 25 
super-pixel array. Clearly, at this distance from the centre of the galaxy, the density is approximately uniform. Fitting a smooth two- 
dimensional model to the data (panel 'b' — fitting details arc given in the text) and subtracting shows the small-scale variations apparent 
in panel 'c'; the significance of these variations is displayed on panel 'd'. 
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Carlo simulations. Artificial stars were added to the frames, 
which were then reprocessed in the same way as the orig- 
inal data, locating the 50% completeness level of the data 
at I ~ 28.8 and V ~ 29.7. The distribution of photomet- 
ric differences between input and recovered magnitudes as a 
function of magnitude shows that the photometry is excel- 
lent with virtually no bias down to very faint magnitudes; 
for instance, at I = 27 the uncertainty is 07 ~ 0.1. We 
have fitted the R05 uncertainty functions with the following 
relations: 

a v (V) = 2.883 x 10~ 9 exp(0.592 , l/ + 0.672) , (1) 
07 (I) = 3.029 x 10~ 7 exp(0.5307 - 1.548) , (2) 

which give a good approximation to the uncertainty esti- 
mates for V < 29 and I < 28. In the analysis below we will 
assume that the I and V uncertainties are uncorrelated (i.e. 
with no colour terms, which is a reasonable assumption), so 
we adopt the above fitting functions to assign uncertain- 
ties to the individual stars in the survey. To investigate 
the properties of the stellar spatial distribution, we select 
only stars with excellent photometry, i.e., (TF6O6W < 0.2 and 
0"F8i4w < 0.2. We further limited the selection to stars with 
—0.5 < (V — I)o < 4.0, yielding a catalogue of approxima- 
tively 43 400 stars. 

As well as encompassing stars in the outskirts of 
NGC 5128, the deep ACS observations also intersect the 
foreground Galaxy. To estimate the foregrou nd counts, we 
used the Besancon Galactic population model (jRobin et al. I 
l2003h to estimate the foreground contamination. The model 
predicts 70-80 stars in the range of I = 21 — 28 over a 
field-of-view comparable to that of ACS in the direction of 
NGC 5128. The right panel of figure [2] shows a typical pre- 
dicted CMD of the Galaxy foreground stars. The model pre- 
dicts about 15-20 Galaxy foreground stars with colours and 
magnitudes similar to those of red giant branch stars at the 
distance of NGC 5128. This constitutes much less than one 
per cent of the total number of stellar objects classified as 
red giant stars in our photometric catalogue, indicating that 
the foreground contamination is not an issue for the analysis 
presented below. 

Large background galaxies, as well as bright foreground 
stars cause holes in the stellar distribution map. To cor- 
rect for this, and for the gap between the two CCDs in the 
camera, when studying the two-dimensional distribution of 
stellar density and metallicity, a mask was constructed by 
choosing suitably large elliptical areas around these prob- 
lematic regions; in all, 3% of the image was discarded in 
this manner. Only those stars outside of the masked regions 
were kept in the final catalogue. 

It is well known that the photometric properties of old 
red giant branch (RGB) stars are linked to their metallic- 
ities, offering therefore a way to determine the metallicity 
distribution of old stellar populations by comparing their 
colour-magnitude diagrams to RGB tracks. The metallicity 
of stars on the red clump cannot be estimated however by 
comparing their photometry to RGB tracks; we therefore im- 
posed a faint cut-off at Mi < —0.75 to the catalogue selected 
above in order to retain only bona-fide RGB stars. This limit 
corresponds to I a < 27.17, or / < 27.38, where the /-band 
measurement uncertainty is 07 = 0.13. An additional qual- 
ity cut of av = 0.13 was imposed to ensure that the V-band 
measurements were also of good quality, yielding a catalogue 



of about 13 300 RGB stars. It is worth mentioning that these 
magnitude cuts do not exclude metal-rich RGB stars. To 
estimate the metallicities of the selected RGB stars, w e re- 
peat the procedure explained in lMouhcine et all J2007|). in - 
terpolating between the models of VandenBerg et al. 1 20061 ) 
of a-enhanced, i.e., [a/Fe]=0.3, RGB stars of mass 0.8Mq. 
The models span -2.014 < [M/H] < -0.097 approximately 
in steps of 0.1 dex, and are complemented at the metal rich 
end by two models with [M/H] = 0.0 and [M/H] = +0.4. 
Stars outside of the colour-magnitude range covered by these 
RGB tracks were flagged and not used in subsequent analy- 
sis (i.e., no extrapolation beyond the validity of the models 
was attempted). At the distance of NGC 5128, the asymp- 
totic giant branch (AGB) bump is found to be located at 
/ ~ 26.8 with a dispersion of 0.12 mag (see R05 for more 
details). The sample of selected bona-fide RGB stars there- 
fore contains AGB bump stars. To examine the effects of this 
on the estimates of the chemical properties of the popula- 
tion of RGB stars, we constructed two samples by selecting 
stars brighter than I — 27.38 and I = 26.6 respectively. The 
metallicity distribution functions of both samples are found 
to be almost identical, i.e., the overall shape is unchanged 
and the changes in the median and /or average metallicities 
are much smaller than the typical estimate of the metallicity 
error. This indicates that the presence of AGB bump stars 
in the selected sample of RGB stars does not bias the results 
discussed below. 



3 ANALYSIS 

3.1 Spatial sub-structure 

We first perform an analysis of the number density variations 
over the area of the ACS survey. Panel 'a' of Fig. shows 
the star-count density map in the X, Y CCD coordinate sys- 
tem chosen by R05. We opted to bin the stars in a grid of 
25 x 25 super-pixel, to contain ~ 100 sources (to ensure a 
signal-to- noise ratio of ~ 10). The star number density has 
been corrected for the mask, and only super-pixels with a 
usable surface area larger than 50% were kept. The map ap- 
pears approximately homogenous, with no obvious gradient. 
This is not very surprising since the observed field covers a 
small physical size at the distance of NGC 5128 compared 
to its projected radial distance. Since we do not have an a- 
priori model of the large-scale stellar distribution over the 
targeted field, we fit the observed stellar distribution with 
a two-dimensional Legendre polynomial function (shown in 
panel 'b'), with up to quadratic terms in X and Y including 
cross-terms (i.e. 6 parameters), which is sufficient to provide 
a smooth empirical template to compare the data to. The 
residuals between the data and the best smooth model are 
shown in panel 'c', along with their statistical significance in 
panel 'd', where we have assumed only Poisson noise uncer- 
tainties (the effect of the mask is properly accounted for). 
There appears to be no coherent large-scale structure in the 
residuals that might indicate the presence of, for instance, 
streams or shells of stellar debris. 

How likely are these small-scale variations about 
the model of the stellar spatial distribution? A way 
to quantify the significance of the presence of sub- 
structure in a survey is to estimate the rms devi- 
ation of the data around a smooth model, taking 
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into account th e expected Poisson noise in the model. 
iBell et all (|2008h proposed the following statistic: <r/total = 
yjk Ei(A - Mi) 2 - I YJM> - E» A, where n 

is the number of discrete bins the sample is divided into, 
Di are the individual observed counts in those bins, Mi are 
the model values in the pixels, and M[ is a Poisson real- 
ization of the model with mean Mi . Thus the numerator is 
the pixel scatter of the data around the model minus the 
expected scatter in the model, while the denominator was 
chosen to be the total of the Poisson scatter. This statistic 
is convenient because it is independent of the choice of bin- 
ning scale (provided that the substructures are well sampled 
by the chosen binning scale), and the number of sources in 
the survey, and furthermore the Poisson noise contribution 
is removed. 

For a bin size ranging from 100 pc to 350 pc, i.e., grids 
of 42 x 42 to 12 x 12 respectively, the fractional rms deviation 
a/total varies from 0.04 to 0.02 with probabilities of finding 
those values of this statistic by chance ranging from 12% 
to 22%, assuming the smooth model shown in panel 'b' of 
Fig-El This clearly indicates that the small-scale deviations 
of the data around the model shown in Fig. [3] are consistent 
with what would be expected from Poisson-noise deviations 
from a smooth model. The stellar populations within the 
surveyed field are then most likely smoothly distributed at 
physical scales larger than ~ 100 pc. 

This contrasts with what was reported for the two large 
late-type galaxies where the properties o f small-sca l e sub - 
structures have been investigated so far. IBell et all (|2008l ) 
found that the fractional rms deviation exceeds er/total = 
0.4 in the Galactic halo region p robed by th e SPS S. For the 
Milky Way analogue NGC 891, llbata et all (|2009l ') have es- 
timated a/total = 0.14 for the extra-planar regions within 
~ lOkpc from the galactic plane, with an extremely small 
probability of finding this value by chance. This clearly in- 
dicates that the spatial distribution of halo stars in both 
large late-type spiral galaxies is most likely dominated by 
significantly lumpy components. 

The catalogue used in the present analysis integrates 
stellar populations along the line of sight through the en- 
tire halo of the galaxy. The contrast introduced by any sub- 
structure will be then reduced, especially in the presence 
of multiple sub-structures along the probed line of sight. It 
should then be appreciated that the low fractional rms de- 
viation found here does not necessarily imply that localized 
spatial over-densities are absent in the halo of NGC 5128. 
Furthermore, the present survey samples a halo region muc h 
smaller than that studied in NGC 891 bv llbata et ail (2009), 
and a tiny fraction of the halo when compared to the vast 
region of the Galacti c halo sampled by the SDSS data ex- 
amined in IBell et al.l (|2008h . Wide-field observations cover- 
ing a representative volume of the outskirts of the galaxy 
are needed to establish (or to rule out) firmly the pre- 
dominance of a smooth stellar component at large radii in 
NGC 5128, and to quantify comprehensively the amount of 
sub-structures in those regions. 



3.2 Metallicity sub-structures 

In this section we investigate the properties of the two- 
dimensional distribution of stellar metallicities over the sur- 



veyed region. Given the smaller size of the sample of stars for 
which we can measure photometric metallicities, we decided 
to distribute stars in a 12 x 12 bin grid (350 pc x 350 pc) to 
ensure a high signal-to-noise ratio in each super-pixel. The 
star number density has been corrected for the mask, and we 
have discarded the few super- pixels where the us able surface 
area was less than 50%. As in llbata et al.l (|2009h . we use the 
median metallicity as a statistic to investigate the spatial 
variations of stellar populations in the outer regions of the 
galaxy. This choice presents the advantages of being insensi- 
tive to e.g. variations of the population spatial density, holes 
due to bright stars, etc. 

The map of median metallicities estimated in the super- 
pixels is shown in the upper panel of Fig. 3] The lower panel 
of the figure shows the map of the uncertainties of median 
metallicities. The errors on the median metallicity in each 
super-pixel is estimated by simulating the effects of both 
the photometric uncertainties on the metallicity estimate 
for each star belonging to the super-pixel and the sampling 
uncertain ty due to the lim ited number of stars per super- 
pixel ( see llbata et alj|2009l . for a detailed description). The 
map shows that the combined random errors on the median 
metallicities in the super-pixels are of the order of 0.02 — 

0. 04 dex. 

The first clear feature of the median metallicity map 
is the presence of a large scale gradient over the field, with 
the average median metallicity decreasing when moving di- 
agonally with increasing Y (i.e., moving from the West to 
the East). A second striking feature of the median metal- 
licity map is that the pixel-to-pixel variations are signifi- 
cantly larger than the random errors. Note that while we 
notice the presence of numerous super-pixels dominated 
by stellar populations with lower median metallicities, i.e., 
[M/H] < —0.7, than that of the overall stellar population, 

1. e., [M/H] ~ —0.6, it appears however that metal-rich, i.e., 
[M/H] > —0.5, super-pixels are rare. How likely are those 
small-scale variations about a large-scale smooth model of 
the median metallicity distribution? 

To represent the large-scale spatial distribution of me- 
dian metallicities, we fitted the observed distribution with a 
two-dimensional Legendre polynomial (up to order 3 in both 
spatial directions, i.e., 10 parameters). The fit is shown in 
the top panel of Fig.0 and exhibits clearly a large-scale me- 
dian metallicity gradient. The middle panel of Fig. [5] shows 
the difference between the map of median metallicities and 
the best smooth model, while the lower panel shows the map 
of the significance levels of these differences. The lower panel 
of the figure shows clearly the presence of localized varia- 
tions at the 3-5 a level in the median metallicity map, occur- 
ring on the scale of a super-pixel or a few super-pixels, and 
distributed over the field. When comparing the significance 
map of the deviations of the median metallicity map about 
the best smooth model to the maps of the spatial density 
distribution and the deviation about the best smooth model 
(see Fig. [3}, we notice that the small-scale sub-structures 
in the median metallicity map are not systematically associ- 
ated to small-scale stellar density sub-structures. As already 
discussed in detail in our study of NGC 891, the number 
density of stars required to cause statistically significant de- 
viations about a smooth median metallicity map is rather 
small, causing only a modest enhancement in the stellar 
distribution that is much harder to detect from its density 
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contrast. The stellar number density needed to cause a de- 
viant super-pixel in the median metallicity map depends on 
the shape of the metallicity distribution functions of both 
the smoothly-distributed population and the super-imposed 
sub-structures, but in broad terms an additional stellar pop- 
ulation numbering 10—15 per cent of the smoot h population 
could induce a ~ 0.1 dex metallicity variation (jlbata et al.l 
2009). While this mass fraction of the super-imposed popu- 
lation is of the order of the change in density over a typical 
super-pixel, the induced change in the median metallicity of 
the super-pixels is many times larger than the typical er- 
rors on the estimate of the median metallicity (see the lower 
panel of Fig-[4j . It is then not surprising that sub-structures 
detected in the median metallicity map are not inducing 
significant deviations in the number density map. 

The distribution function of the statistical significance 
of the variations about the best smooth model of the me- 
dian metallicities map is shown in Fig. [6] The dashed line 
shows a fitted Gaussian of dispersion equal to 1.58. The 
probability of finding such a large dispersion by chance (we 
expect unit dispersion if the uncertainties have been cor- 
rectly estimated), given the large number of bins, is smaller 
than 0.1%. This strongly indicates that the variations about 
the best smooth model of the spatial distribution of me- 
dian metallicities are significantly larger than expected for 
a population where the median metallicity varies spatially 
in a smooth manner. 

Perhaps the simplest explanation of these large devia- 
tions is that we have underestimated the uncertainties in 
the median metallicity by 30 — 50%. This would require 
that the photometric errors reported by R05 to be underesti- 
mated by a similar factor. In our study of NGC 891 we were 
able to verify our uncertainty estimates by comparing mea- 
surements in overlapping fields. Unfortunately, we do not 
have the necessary complementary data to perform a simi- 
lar check here. Although we cannot dismiss this possibility, 
experience from our NGC 891 study (and other projects) 
suggests that artificial star simulations give reliable uncer- 
tainty estimates. 

The large pixel-to-pixel variations reported in the me- 
dian metallicity map are unlikely to result from errors in 
the zero-point over the ACS field. To induce the observed 
small-scale metallicity changes, large variations in the zero 
point (of the order of 0.1 mag or more) are required, which 
is extremely unlikely for the very stable instrument that the 
ACS is. 

Could the large-scale metallicity gradient and the small- 
scale chemical substructures we have uncovered be due to 
systematic biases? To convert the RGB photometric proper- 
ties into metallicities, we have had to assume that the stellar 
populations dominating the surveyed field are uniformly old. 
Due to the so-called age-metallicity degeneracy, a red RGB 
star can be either more chemically evolved than a bluer one 
or older. In the interpolation procedure, we have assumed 
a single old age for all the stars, which may introduce a 
bias of up to 0.1 de x, if the age is ~ 4 Gyr younger (e.g. 
iReikuba et"ai~ll2005l ). 

The AGB stars are known to be efficient trac- 
ers of intermediate age ( 1-6 Gyr) stellar popul ations 
l|Mouhcine fc Lancon ll2002l ; iMarigo fc Girardi \\200H ). The 
abundance of AGB stars normalized to RGB stars, which 
are tracing predominantly older populations, is a good in- 



dicator of a stellar population age (e.g. iFrogel et al~lll990l ). 
The map of the ratio of AGB stars, selected as stars with 
22.0 < Io < 23.5 and 1.5 < (V — I) < 3.0, and all other 
quality criteria identical to the RGB sample summarized 
in Sj2] to stars in the brightest one magnitude of the RGB 
is shown in Fig. [7] The spatial distribution of the AGB to 
RGB ratio is fairly smooth, with no apparent spatial gra- 
dient. This suggests the absence of an AGB gradient and 
then most likely the absence of an age gradient across the 
observed field. Furthermore, no enhancements of AGB stars 
relative to RGB stars are associated to any of the super- 
pixels that are deviants from a smooth model of the median 
metallicity map by more than 3a. Interestingly, the typical 
AGB to RGB r atio, i.e.. < 0.2 . is consistent with an old stel- 
lar population (Mould 2005). The pixel-to-pixel variations 
of the AGB-to-RGB ratio are significantly smaller than what 
is required for the pixel-to-pixel variations of the median 
metallicities to be due only to age differences: to account 
for the > 0.15 dex variations in median metallicity about 
the best smooth model for the most deviant super-pixels 
(i.e., ±3cr or more) as a result of age differences, the stellar 
populations of those super-pixels have to be dominated by 
< 6 Gyr stars, which is not suppor ted by the morpholog y of 
the colour-magnitude diagram (see IReikuba et al. I[2005l . for 
a detailed discussion). 

The localized median metallicity variations are also un- 
likely to result from the effects of small-scale structures 
in the foreground Galactic inter-stellar medium. Following 
llbata et al.l l|2009l ). by using the reasonable assumption that 
reddening variations scale approximately linearly with mean 
reddening, we expect less than a couple hundredths of a 
magnitude scatter at most in colour excess from foreground 
dust. For the small-scale metallicity variations to be due 
solely to variation in foreground extinction, the required off- 
set in colour excess are either much larger than the expected 
scatter, or, more problematically, negative in the case for de- 
viant metal-poor super-pixels. Small angular scale variations 
in foreground dust are therefore unlikely to contribute sig- 
nificantly to the observed small-scale variations in median 
metallicity. 

We conclude that there is no evidence for systematic 
biases affecting the measure of stellar metallicities. So, un- 
der the assumption that the photometric uncertainties have 
been well estimated, we conclude that the statistically sig- 
nificant small-scale variations in the median metallicities of 
the stellar populations are present in the targeted halo field. 
To highlight the differences between thte detected small- 
scale chemical inhomogeneities, Fig. [5] shows the CMDs of 
all stellar objects populating a metal-poor super-pixel, i.e., 
with a median metallicity of [M/H]=-0.77, and a metal-rich 
super-pixel, i.e., with a median metallicity of [M/H]=-0.47, 
located at (X/kpc=1.75, Y/kpc=3.00) and (X/kpc=2.40, 
Y/kpc=0.48) respectively. Fiducials of RGB sequences for 
Galactic globular clusters spanning a wide range of metal- 
licities and shifted to the distance of NGC 5128 are over- 
plotted. The stellar populations occupying both super-pixels 
appear to be different. A significant fraction of bright RGB 
stars populating the metal-poor super-pixel exhibit photo- 
metric properties consistent with metallicities in the range 
of [Fe/H]=-2.2 and -1.3. Such metal-poor stars are absent in 
the metal-rich super-pixel. Furthermore, the reddest stars 
in the metal-rich super-pixel extend clearly to higher metal- 
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Figure 4. Panel 'a' displays the median metallicity calculated 
in each of the spatial super-pixels. A slight gradient is apparent, 
with the top left-hand corner being more metal-poor than the 
bottom right. The corresponding uncertainties, calculated from 
Monte-Carlo simulations accounting for both the photon noise 
and the sampling errors, are shown on the panel 'b'. 

licities than those occupying the metal-poor super-pixel. 
The stellar population differences between the localized sub- 
structures in the median metallicity spatial distribution 
would suggest that the survey has revealed the presence of 
inhomogeneities that have not yet been fully blended into 
the smooth halo. 



4 DISCUSSION 

The metallicity distribution function of the stellar popu- 
lations in the outer regions of NGC 5128 has been deter- 
mined over an exte n ded range of radial distances, 10-40 kpc 
(IHarris et al ] 1 19991 ; IHarris fc Harrisl |2002| ; iReikuba et al. I 
l2005h . The overall shape, mean metallicity, metallicity dis- 
persion, fracti on of metal-poor star s, are found to be strik- 
ingly similar (|Reikuba et al. 1 120051 ) over these fields. This 
indicates that the metallicity gradient present within the 
field at ~ 40 kpc (see Fig. [5]) could not be associated to a 
large-scale metallicity gradient, as it would lead to metal- 
licities significantly higher than observed if extrapolated in- 
ward. The most likely explanation is that the median stellar 
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Figure 5. A two-dimensional polynomial model fit to the median 
metallicity map (panel 'a' of Fig.[4j is shown in panel 'a'. The 
residuals between the data and this model are displayed in panel 
'b'. Using the uncertainty estimates on the metallicity metallicity 
from panel 'b' of Fig. [4] the statistical significance levels of the 
metallicity variations are shown in panel 'c'. 



metallicity gradient present within the ACS halo field is re- 
lated to chemical inhomogeneities present locally. 

Using the predictions of a A-cold dark matter semi- 
analytical galaxy formation model, and by summing the 
total star formation occurring over all the progenitors, 
iBeaslev et al. I (|2003l ) computed the metallicity distribution 
function of spheroid stars of early-type galaxies comparable 
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Figure 6. The distribution of statistical significance of the de- 
viations from the best smooth model of the median metallicity 
map shown in the top panel of Fig. [5] The Gaussian model fit 
superimposed on this distribution has a dispersion of 1.58. The 
distribution of significance levels has an extremely small proba- 
bility of being drawn by chance from a Gaussian distribution of 
unit dispersion. 
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Figure 7. The spatial distribution of the ratio of AGB stars to 
RGB stars over the ACS field. 



in luminosity and environment to NGC 5128. Despite some 
differences, the predicted metallicity distribution function 
resembles qualitatively the observed distribution function 
of stellar metallicities in the outer regions of NGC 5128, 
i.e., mean metallicity, m etal-rich end cut-off, and extent of 
the metal-poor tail fsee iReikuba et al. 1 120051 ') . This broad 
agreement suggests that the stellar populations in the outer 
regions of NGC 5128 could have been assembled in a hier- 
archical fashion. If this is indeed the case, one can expect to 
find the presence of chemical sub-structures caused by the 
remnants of the disruption of the predicted large popula- 
tions of sub-halos orbiting around massive galaxies. 

The localised chemical inhomogeneities could be in prin- 
ciple potential candidate satellites or remnants of past ac- 



cretion events. Upon inspecting the stacked ACS image 
visually, none of the strongly deviant super-pixels shows 
any plausible stellar concentration, although clear RGB se- 
quences are detected for all of them, apart from the super- 
pixel at (X ~ 1.8 kpc, Y ~ 3kpc) which contains the faint, 
old, and metal- p oor ex tended globular cluster reported in 
iMouhcine etafl (|2010bl) . The extended nature and the low 
luminosity of that globular cluster become very important 
when considered in light of the emerging consensus that ex- 
tended chistere_are fOT^ mass galax- 
ies (e.g |Elmegreen~ll2008l ; lDa Costa et al. Il2009t ). In M31, all 
but one extended globular cluster in the outer halo, i.e., at 
projected distances larger than 30 kpc, are found to be pro- 
jected onto s tellar substructures o r are members of a cluster 
overdensity (|Mackev et al. I [20101 ). supporting the scenario 
where this class of globular clusters originate in low mass 
systems. 

The properties (luminosity, surface brightness, size, cen- 
tral concentration) of the globular cluster found in the halo 
field resemble in every way one those of the faint and ex- 
tende d globular clusters pop ulating the Milky Way outer 
halo (|Mouhcine et alj l2010bj ) . This suggests that the ex- 
tended cluster in the targeted field could have been formed in 
a similar environment as its Galactic counterparts, and ex- 
perienced potentially a similar dynamical history. Strong ev- 
idence is accumulating that extended and faint globular clus- 
ters in the outer Galactic halo, i.e., the so-called young halo 
globular clusters, were ass ociated with dwarf galaxies that 
have since disrupted (e . g. Zinnl[l9 93; Mackcy & Gi lmorel 
|2004 iForbes fc Bridges! l2010t ). iMackev et al. I l|2010h have 
argued that the vast majority of globular clusters in the 
outer halo of M31 are most likely associated with tidal de- 
bris features observed in those regions, thus strongly sup- 
porting the scenario where the outer globular cluster system 
of M31 has been built up via the accretion of satellite host 
galaxies. It is then tempting to suggest that the extended 
faint cluster in the surveyed halo field in NGC 5128 could 
have originated in a now disrupted dwarf galaxy. Interest- 
ingly, the median metallicity of the diffuse stellar population 
within the super-pixel containing the extended faint clus- 
ter is comparable to those measured for metal-poor super- 
pixels deviants at the 3 — 4 a levels in the observed halo field, 
i.e., [M/H] ~ —0.75, suggesting that the stellar populations 
dominating the metal-poor super-pixels could be sharing the 
same origin. Put together, these properties provide circum- 
stantial evidence that the chemical inhomogeneities detected 
in the observed field could be related to the disruption of low 
mass satellites. 

The restricted spatial extent of the present survey lim- 
its however our ability to investigate the large-scale distribu- 
tion of the chemical sub-structures, to quantify the degree of 
spatial sub-structures in the halo of the galaxy, to search for 
signs of the disruption of the hypothetical progenitor (s) of 
the detected chemical sub-structures. Extending the survey 
of the outer regions of NGC 5128 should be undertaken to 
tackle comprehensively these questions. When compared to 
the properties of sub-structures in the outskirts of late-type 
galaxies, it will be possible to constrain the assembly histo- 
ries of galactic halos as a function of galaxy e.g. morphology, 
type, and mass. 
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Figure 8. CMDs of stars populating two super-pixels with their respective positions indicated in each panel. The left panel shows 
the CMD of the most metal-poor super-pixel in the map shown in Fig. [4] with a median metallicity of [M/H]=-0.77, while the right 
panel show the CMD of a typical metal-rich super-pixel with a median metallicity of [M/H]=-0.41. Note that the median metallicity 
for the whole field is [M/H1=-0 .53. Shown also as solid lines are the red giant sequences of the indicated Galactic globular clusters from 
|Pa Costa fc Armandroff assuming an intrinsic distance modulus for NGC 5128 of 27.9. It is clear that the metal-poor side of the 

red giant branch is much more populated in the case of the metal-poor spatial super-pixel than for the metal-rich one. 



5 SUMMARY & CONCLUSIONS 

To make the first attempt to quantify the degree of substruc- 
ture in the outer regions of elliptical galaxies, we have under- 
taken a structural and chemical analysis of the stellar con- 
tent of a 3.8 kpc x 3.8 kpc halo field situated ~ 38kpc from 
the centre of the nearest giant elliptical galaxy NGC 5128. 
The very deep colour-magnitude diagram of the stellar pop- 
ulations used in the analysis reaches the red clump. 

Comparing the observed stellar number density map of 
halo stars to smooth models, we find no statistically signifi- 
cant evidence for the presence of spatial density variations. 
Deviations from the smooth models, quantified using the 
rms deviation of the data around the best fit model are con- 
sistent with what would be expected from Poisson-noise de- 
viations in a smooth model. The measured statistics appears 
to suggest the presence of a smooth stellar halo component 
at spatial scales larger than ~ 100 pc in the region probed 
by the present survey. 

Estimating metallicities of RGB stars from their op- 
tical photometry, we investigated the spatial distribution 
of stellar metallicities. On top of a smoothly distributed 
large-scale map, we detect the presence of small-scale sub- 
structures in the median metallicity across the surveyed 
field. Although the absolute median metallicity variations 
are relatively modest, i.e., ranging from ~ 0.08 to ~ 0.2 dex, 
they are nevertheless statistically significant. Those small- 
scale chemical sub-structures do not appear to be associated 
with any statistically significant stellar spatial density vari- 
ations. We argue that these subtle small-scale metallicity 
variations are related to accretions of low mass systems. 
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